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Abstract 

Motivated by the absence, so far, of any direct signal of conventional low-energy super- 
symmetry, we explore the consequences of making the lightest Higgs boson in supersymme- 
try relatively heavy, up to about 300 GeV, in the most straightforward way, i.e. via the 
introduction of a chiral singlet S with a superpotential interaction with the Higgs doublets, 
XSH1H2. The coupling A dominates over all the other couplings and, to maintain the suc- 
cessful perturbative analysis of the ElectroWeak Precision Tests, is only restricted to remain 
perturbative up to about 10 TeV. The general features of this "ASUSY" framework, which 
deviates significantly from the MSSM or the standard NMSSM, are analyzed in different ar- 
eas: ElectroWeak Precision Tests, Dark Matter, naturalness bounds on superparticle masses, 
and LHC signals. There is a rich Higgs/Higgsino sector in the (200 - 700) GeV mass region, 
which may include LSP Higgsino dark matter. All other superpartners, apart from the top 
squarks, may naturally be heavier than 1-2 TeV. This picture can be made consistent with 
gauge coupling unification. 



1 Introduction 



The naturalness problem of the Fermi scale in the Standard Model (SM) amounts to understand 
the lightness of the Higgs boson relative to any mass scale in the theory, whatever it may be, that 
completes the SM itself at high energies, or in the Ultra- Violet (UV). If the relative Higgs lightness 
is not accidental, this gives the best hope for seeing new physics at the LHC, hence the focus on 
the naturalness problem. As is well known, the constraints on a positive solution to this problem 
largely come from the success of the SM in accounting for the results of the Electro Weak Precision 
Tests (EWPT). In fact, the difficulty to solve it is exacerbated by the apparent lightness of the 
Higgs boson, as indirectly implied by the EWPT [1]. If, for one reason or another, we were misled 
in interpreting the data of the EWPT in terms of a light Higgs mass and the Higgs boson were 
heavier, the upper bound on the naturalness cutoff of the theory of the electroweak interactions 
might be relaxed by a non-negligible amount [2, 3]. 

This simple observation has a counterpart in supersymmetry, which is in many respects the 
most natural UV completion of the Standard Model (SM). In supersymmetric extensions of the 
SM, like the Minimal Supersymmetric Standard Model (MSSM), there is no problem in letting the 
cutoff become arbitrarily high, while the lightest Higgs boson is actually predicted to be light. So 
light in fact that to make its mass consistent with the experimental lower bound (about 114 GeV 
in most of the parameter space [4]) requires a large radiative correction due to a heavy stop, which 
in turn brings back at least a few percent fine-tuning in the Z boson mass. Accommodating a 
heavier Higgs boson would therefore ease the naturalness problem in supersymmetry too, while 
allowing heavier superpartners. This has motivated several attempts in the literature to increase 
the lightest Higgs boson mass [5 - 16]. 

Here we explore the possibility of making the lightest Higgs boson in supersymmetry relatively 
heavy in the most straightforward way, i.e. via the introduction of a chiral singlet 5* with a large 
superpotential interaction XSHiH 2 , where Hi are the two Higgs doublet multiplets, leading to a 
large quartic interaction of the Higgs doublets. Unlike the usual treatment of the Next to Minimal 
Supersymmetric Standard Model (NMSSM) [17], A is not perturbative up to the unification scale. 
We insist only that consistency with the EWPT can be analyzed in a fully perturbative manner, 
which we ensure by requiring that A remains perturbative up to pa 10 TeV. At this scale a change 
of regime of the theory should intervene, which we leave unspecified. It is an open question 
whether this will allow a successful unification of the gauge couplings. A positive existence proof 
of this is given in Refs. [10, 11, 12]. 

With perturbativity preserved up to ~ 10 TeV, contributions to the EWPT from unknown 
UV physics can be sufficiently small. There may, however, still be calculable "infrared" contribu- 
tions which make the theory incompatible with the EWPT. In particular, one might think that 
contributions from a lightest Higgs boson of mass (200 ~ 300) GeV, as well as those from the other 
Higgs sector particles which grows like the fourth power of A, are dangerous for the EWPT. Quite 
on the contrary, our general analysis of the contributions from the Higgs boson sector and, even 
more importantly from the Higgsino sector, shows that the theory is perfectly consistent with the 
EWPT in a significant range of the parameter space. Deviations from the "EWPT ellipse" due to 
a relatively heavy lightest Higgs boson are entirely canceled by the contributions from the other 



1 



Higgs bosons and the Higgsinos, in almost all the parameter space in which electroweak symmetry 
breaking occurs naturally. 

In Section 2 we define our theory, minimize the scalar potential for large A, and give both 
the scalar and fermion spectrum. Many significant features of the theory are independent of the 
specific forms for the superpotential and soft operators, but depend critically on A being large. 
Hence we call this framework "ASUSY." In Section 3 we consider the EWPT from both scalar 
and fermion sectors. In Section 4 we consider the candidate for Dark Matter (DM) in the special 
case of heavy weak gauginos, as allowed by naturalness. In Section 5 we obtain the naturalness 
bounds on the various particles in the limiting case in which the scalar S is taken heavy. The 
characteristic manifestations of the model at the LHC are described in Section 6. Summary and 
Conclusions are given in Section 7. 



2 The ASUSY Model 

We consider the most general supersymmetric theory with SU(3)c X SU{2) L x U(1)y gauge 
invariance that contains a singlet chiral superfield S in addition to the fields of the MSSM. The 
theory possesses the superpotential interaction 

W X = XSH 1 H 2 , (1) 

where Hi and H 2 are the usual Higgs doublets, coupled respectively to the down-type and up- 
type quarks. We increase the mass of the lightest Higgs boson, and therefore the naturalness 
of the theory, by taking A large. We call this the ASUSY model, to stress the key role of the 
large coupling A. Such a large coupling may arise from compositeness of (some of) the Higgs 
states [10, 11, 12, 15]. 

How large can A be? To keep perturbativity up to at least 10 TeV, we require A 2 (A = 10 TeV) 
be less than about 4tt. Given the Renormalization Group Equation (RGE) satisfied by A (t = In E) 

A^ m 
dt ~ 2tt 2 ' [Z) 

this implies that, at the low-energy scale of about 500 GeV, A is less than about 2. Throughout 
the paper, unless differently stated, we take A(500 GeV) = 2. With this value of A the Landau 
pole is typically a few tens of TeV. 

The superpotential of our model takes the general form 

W = t i(S)H 1 H 2 + f(S), (3) 

and, neglecting the gauge terms, the scalar potential can be written in the form 

V = rf(S)\Hi\ 2 + fi 2 2 (S)\H 2 \ 2 - (4(S)HiH 2 + h.c.) + X 2 \HiH 2 \ 2 + V(S). (4) 

The mass parameters /if (s), /^(s), /^(s), /x(s) and 

M^iM, ,5) 
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calculated at the background expectation value s of the field S 1 , characterize many of the properties 
of the theory. We generally denote them as /if , fi 2 , fj%, fi, and M, by leaving understood their 
argument s. For simplicity we assume CP invariance of V and W. This formulation of the ASUSY 
theory is convenient because many equations take a similar form to the familiar MSSM case. 
The stability of the potential (4) requires 

pi, ill > |/i| 2 , (6) 

for all the values of s, whereas the condition for electroweak symmetry breaking is 

\l4\ > (JiifJto, (7) 

at the minimum of the effective potential for s, obtained by replacing Hi and H 2 by their s- 
dependent expectation values. Under the condition (6) electromagnetism is unbroken. To be 
more precise, once the gauge contribution to the potential is introduced, the potential is always 
stable, and the condition (6), with an additional small term on the right-hand-side, becomes the 
condition for unbroken electromagnetism. For our purposes, this difference has negligible impact. 
In general we neglect in this section all the relatively small effects of the SU(2)l x U(1)y D-term 
contributions. We shall also neglect the possibility of spontaneous CP violation, i.e. we shall 
consider all the mass parameters real. 1 

By minimizing the potential it is straightforward to find 

tan/3 = ^ = ^, (8) 

Vl 111 

\2 2 ^3 2 2 (r\\ 

where A = dfi(S)/dS, and V\, v 2 , v = (v \ + v^) 1 ^ 2 are the usual Vacuum Expectation Values 
(VEVs) of the Higgs fields. One also finds that the mass of the charged Higgs bosons, H ± , is 

m 2 H± = n\ + fj%, (10) 

where [i\ and fi 2 , should be evaluated at the minimum of s. Unlike the mass of H ± , the masses 
of the neutral scalars h, H and A can be expressed in terms of f/,i, (/,% and only if their mixing 
with the scalar S can be neglected. In this case, the mass of the pseudoscalar Higgs state is given 
by 

while the two CP-even Higgs states hi = y/2(ReHf — Vi) mix with the following mass matrix: 



m\ sin 2 /? (X 2 v 2 - hm 2 A ) sin 2(3 



2 

(X 2 v 2 — \m 2 A ) sin 2(3 m 2 A cos 2 (3 



(12) 



1 A physical phase in the Higgsino mass matrix (sec below) could have relevant consequences. 
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where ii\ in Eq. (11) should be evaluated at the minimum of s. The masses and composition of 
the mass eigenstates H and h are given by 



H,h 



and 



respectively, where 



m\±X), X 2 = m\- AX 2 v 2 m 2 H± sin 2 2(3, (13) 



H = cos a hi + sin a h 2 , (14) 
h = — sin a hi + cos a hi, (15) 



m 2 A cos2(3 + X 

tana= (AV-< ± )sin2/r (16) 



We find that the following inequalities hold: 

rrih < sin (3 Xv, mn± < ran < tua- (17) 

In the fermion sector, the SU(2)l x f/(l)y gauginos are mixed with the other fermions only 
by the relatively small gauge terms. Since there is also no significant naturalness upper bound on 
their masses, it is relevant to consider the case in which these gauginos are decoupled from the 
other fermions. In this case the charged Higgsino has a mass [i and the mass matrix of the neutral 
fermions 

Ni = ^=(Hi-H 2 ), N 2 = ^=(Hi + H 2 ), N 3 = S, (18) 
is given in terms of \x and M by 

f H A- 



M 



V2 

-A^p 

vA ^7r - a "tt m J 



(19) 



where /x and M are, again, evaluated at the minimum of s. Only the relative sign of \x and M 
is physical, and we fix our convention to be \x > for definiteness. Note that the condition (6) 
implies the inequality 

ji < cos j3 m H ±. (20) 

The iVj's will be related to the mass eigenstates (neutralinos) \% by a non-trivial 3x3 mixing 
matrix V: 

Ni = V ijXj , V T V = l, V T MV = diag(m 1; m 2 ,m 3 ). (21) 
The lightest neutralino x has a mass satisfying 

|m x |</i, (22) 
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and the other two eigenvalues are above \i and below — /x, respectively. Moreover, the determinant 
of the mass matrix, and thus m x , vanishes for a particular value of M: 



m x = for M = sin 2(3. (23) 

Thus there is a region of parameter space with a naturally light neutralino, which is an admixture 
of S and the Higgsinos. In fact, the lightest neutralino is much lighter than the characteristic 
scales for fi and M in a wide region around the point satisfying Eq. (23). 

The mixing matrix and the spectrum are difficult to determine analytically, except in the 
limiting case of tan/3 = 1, which leads to 

mi = /i, (24) 
m 2 , 3 = ^(M - /i±Y), Y 2 = (M + /i) 2 + 4AV. (25) 

In this case Xi = -^i> while X2,3 are mixtures of A^ 2 ,3 with a mixing angle 7, where tan 7 = 
— (M + fi + Y)/2\v. The lightest state is either xi 01 X2 depending on whether M is above or 
below a critical value M c = fi — X 2 v 2 /2fi. 

For sufficiently small tan (3 — 1, one can approximately diagonalize the matrix (19) expanding 
in .M13. The mixing between N\ and will be small everywhere except near M = M c . As 
a result, the coupling of the lightest neutralino x to the Z boson will be generically suppressed. 
This may play an important role in the interpretation of DM, as discussed in Section 4. 

When /ii = /i 2 , giving tan/3 = 1, an SU(2) custodial symmetry survives, and plays an impor- 
tant role. In this limit, in the scalar sector either h or H, depending on the sign of 2m 2 H± — m\, 
behaves like the SM Higgs boson, of mass rrih = \v/\/2, whereas the other scalar with H ± forms a 
degenerate S'f/(2)-triplet. Similarly, in the fermion sector, the chargino and Ni form a degenerate 
triplet of mass fi. 



3 Electro Weak Precision Tests 

We perform the analysis of the EWPT in the usual S-T plane, with the experimental contours 
taken from [18]. With respect to the SM, the ASUSY model has: 1) a modified contribution to 
S and T from the scalar Higgs sector [10, 16]; 2) new and relevant contributions from the stop- 
sbottom and from the Higgsinos (we consider heavy gauginos). The contributions of 1) and 2) 
can be analyzed separately. 

To locate our model in the S-T plane before the stop-sbottom and Higgsino contributions are 
added, we subtract the one-loop SM Higgs contributions 

T mggs (m h ) = -3 [A(m h , m w ) - A(m h , m z )] , (26) 
Smgg S {m h ) = F(m h , m z ) + m 2 z G(m h , m z ), (27) 
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Figure 1: The black (darker) curve shows the SM results with a Higgs mass nth = 100 — 350 GeV 
in 50 GeV increments. The ellipses show the regions of the S-T plane allowed by EWPT at lcr 
and 2a. The red (lighter) curves give the predictions from the Higgs scalar sector in ASUSY, 
as described in the text, with values of tan/3 in the interval tan/3 = 1 ~ 5 as indicated and 
m H± = 350, 500, 700 GeV. 

from the SM values of 5* and T (see Appendix A for the definitions), and then add the one- loop 
2 Higgs-Doublet Model (2HDM) contributions (computed under the assumption of no doublet- 
singlet mixing) [19, 20] 

T 2 rdm = sin 2 (/3 - a) [T HiggfS (m h ) + A(m H ±,m H ) - A(m A , m H )\ 
+ cos 2 (/3 - a) [T mggs (m H ) + A(m H ±,m h ) - A(m A ,m h )] 
+ A{m H ±,m A ) : (28) 

SWdm = sin 2 (/3 - a) [S mggs {m h ) + F(m A , m H )\ 

+ cos 2 (/3 - a) [S mggs (m H ) + F(m A , m h )] 

-F(m H ±,m H ±). (29) 

Here, in Eqs. (26, 27) represents a reference value of the Higgs boson mass in the SM, while 
that in Eqs. (28, 29) is the mass of the lightest CP-even Higgs boson in the ASUSY model. 

In the ASUSY model, the 6 parameters appearing in Eqs. (28, 29), i.e. m^, m#, m A , m H ±, a 
and P, depend only on tan/3, and A, which we set equal to 2. In Fig. 1 we show the result for 
the EWPT for several values of tan (3 and m#± (tan (5 — 1~5 and m H ± = 350, 500, 700 GeV). Two 
features are manifest from this figure: i) the role of the custodial symmetry for tan j3 approaching 
unity, thus suppressing the corrections to T; ii) the fact that the positive T-correction brings most 
of the points of the ASUSY model inside the region preferred by experiments, at least as long as 
tan j3 is not too large. 
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Figure 2: The expected stop-sbottom contribution to the T parameter is constrained by natural- 
ness to lie above the curve. 

The stop-sbottom contributions in the zero left-right mixing limit are given by 

,4 / rnn n„\r\ 2 



T s , sb = 6A(m- tL ,m- bL ) « - m * - 2 « 0.05 (^) , (30) 
L 32ir' l a era v z m~ \ m~ tL J 

s N 1 m, 2 /500 GeV\ 2 , , 

5st- sb = FK i ,m^)-FK- i ,m t -J^- — -f ^-0.003 J , (31) 

where the approximate expressions follow from m 2 ^ — m~ = m 2 in the large m 2 ^ limit. While 
the contribution to S is always negligibly small, this is not the case, as is well known, for the 
contribution to T. Anticipating a tan /^-dependent upper bound on the stop masses from natural- 
ness considerations (see Section 5), we show in Fig. 2 the minimum value of T st _ s b when the stop 
masses are taken at this boundary. 2 This strongly reinforces the conclusion that tan (3 cannot be 
too large. If the Higgs and stop-sbottom sectors were the only contributions to S and T, any 
value of tan (3 above 5 or so would be almost excluded. 

We finally consider the contributions from the Higgsino sector described in the previous Section. 
The couplings to the gauge bosons are 

Ant = 9 fV; (-V-fNx + ®Yl 5 N 2 ) + h.c. 

+ \K (* 7 m * + NiYl 5 N 2 ) + (ty** - N x ^-fN 2 ) , (32) 

where \& represents the chargino. The contributions of Higgsinos to T and S can be written in 



2 Taking into account the left-right mixing may help reduce T s t_ s b, although not dramatically, because the 
corresponding A-term parameter, At, will be subject to the same naturalness bound as wlq and mi R . 
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terms of the mixing matrix V as follows: 



Tffiggsinos =^ a=i (Vi a ) 2 A(/i,m a ) + (V 2a ) 2 A(n, -m a ) 

-\Y! \ {V la V 2b + V lb V 2a ) 2 A{m a ,-m h ), (33) 

X V ^3 ry ~ ~ 

Sffiggsinos =7T / J L (^laV 2 6 + VlfcV^a) T(m a , -™&) ~ -F(/i, A*)) (34) 
Z ' 'a, 0=1 

(see Appendix A for the definitions of the functions A and F). These contributions are shown 
in Fig. 3 in the \x-M plane. Also shown in the plot is the region (shaded) in which the lightest 
neutralino mass is less than half the Z boson mass: m x < mz/2. Even though its coupling to the 
Z boson is suppressed, this mass range is inconsistent with the negative searches from the LEP. In 
view of Figs. 1 and 2, as well as the mostly positive contributions to S and T in the plots of Fig. 3, 
we conclude that most of the \x-M plane is allowed for tan/5 ^ 3, except for a tan /^-dependent 
strip around M = where T becomes too large. 

The dependence of Tkiggsinos on tan (3 can be approximated by 3 

T H ig gsi nos ~ -F(tan/3) T(/i, M), T(t) = (h^j) > ( 35 ) 

where T(/i, M ) is some function of \x and M. To make this scaling evident, the values for the 
contours of T in each panel of Fig. 3 are chosen such that they are proportional to the contours 
of the tan/3 = 1.5 panel, scaled by a factor jF(tan/3)/jF(1.5). Even for tan/5 > 3 there remain 
two strips around M ~ 600 GeV and M ~ —800 GeV where THiggsinos is sufficiently small to be 
allowed (especially combined with the positive S parameter). However, taking into account the 
contribution T st _ s b, which starts being increasingly problematic for tan/5 (4~5), as well as the 
fact the allowed region for \i from naturalness rapidly shrinks for tan/3 ^ (4 ~ 5), we conclude 
that the region 

tan (5 ^ 3, (36) 

is preferred in the ASUSY model. 



4 Dark Matter for Heavy Gauginos 

There are several possibilities in ASUSY for the Lightest Supersymmetric Particle (LSP), corre- 
sponding to differing LHC signals and DM candidates. In this section we study the particular 
case that the LSP is dominantly Higgsino, having only small gaugino components. In the previous 
section we have shown that decoupling the gauginos leads to a large region of parameter space 

3 This relation can be understood easily in two limiting cases: 1) for tan/3 sa 1, in which case we can compute 
T perturbatively expanding in the M13 entry of the Higgsino mass matrix. 2) in the decoupling limit of heavy 
Higgsinos. in which case integrating them out produces a dimension 6 operator \H\D^Hi — li\D ^H^^ ■ However, 
Eq. (35) apparently holds with reasonable accuracy even in the general case of light Higgsinos and/or tan/3 — 1 
being not so small. 
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Figure 3: The contributions from the Higgsinos to the T parameter (red color; contour values not 
framed) and the S parameter (blue color; contour values framed). Vertical dashed lines represent 
an upper limit on \i coming from naturalness. Shaded regions indicate the regions in which the 
lightest neutralino is lighter than rriz/2 ~ 45 GeV. 
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consistent with EWPT, and in this section we find that the presence of the singlet Higgsino allows 
the interesting possibility of Higgsino DM, whereas in the MSSM it is excluded. 

We now compute the thermal relic abundance of the lightest Higgsino x- We use the standard 
formalism from Ref. [21]. The freeze-out point is given in terms of the scaled inverse temperature 
x = m x /T: 

0.038 g x m Pl m x {av rel ) 

x f -In ' W 

9* %f 

where m P \ ~ 1.22 x 10 19 GeV, g x = 2, is the number of SM degrees of freedom relativistic 
at the time of freeze-out {g* = 86.25 for m b <C Tf ^ m w ), and (crv Te i) is the thermal averaged 
annihilation cross section (% + x ~~ ^ a U), which enters the Boltzmann equation for the x number 
density n: 

| + 3ffn = -K)(n 2 -ny. (38) 
At freeze-out, x is nonrelativistic and the cross section can be expanded in v. 

(TV xe \ = a + fefrel, (c^rel) = CL + 65/ X. (39) 

Defining the annihilation integral 

T . , f°° (av rel ) a 3b 

J(x f ) = / '—^dx « - + -j, (40) 



the present-day mass density of x is given by 



_ u2 1.07 x 10 9 GeV" . A1 . 

Q xh 2 * —tj- 2 — — ■ (41) 

gj mpiJ(Xf) 

There are several caveats in the analysis described above. 1) We ignore possible coannihilations. 
This will be valid if all the other neutralinos and charginos have masses at least a fewxTy higher 
than m x . This treatment is indeed justified in the parameter region of interest to us, but may 
not be in general. 2) We do not give a careful treatment of near-threshold situations. Indeed, 
there are important thresholds associated with the W and Z bosons. Near these thresholds, our 
calculation may be subject to larger errors, and a more careful treatment should proceed along 
the lines of Ref. [22]. 

Let us now calculate various annihilation processes in turn. We first calculate x + X ~ * // y i a 
s-channel Z exchange. The x coupling to Z is given by 

~^ Z&f-fx, *x = W lx V 2x . (42) 

The resulting annihilation cross section in the approximation of massless final states is atfV ie i = 

b ff = 7.31/4 9 r7 - o m * 9NO , (43) 
ff x 32ttc4 (4m2 - m|) 2 ' V ' 
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where 7.31 is the sum of g\ + g\ over all SM quarks and leptons except for the top quark. The 
coupling k x is generically suppressed, especially in the lower right corner of the fi-M plane. For 
this reason, large values of Q x h 2 are attained for modest values of m 



v 



Above the W thresholds, we have to take into account the annihilation into W + W~ proceeding 
via s-channel Z exchange and i,ti-channel chargino exchange. The amplitudes for these processes 
are given by 



9 < 



mi 



x [rTik. - k + ) x - r^(k+ + 2k^ + ^(2k + + *_)"] <(fc + K(M, 



and 



9_ 
4 



VjfP 



(k+- pi) 2 



P = Vi 



(44) 



(45) 



respectively, where u\ t 2 are the initial state spinors with momenta pi^, Uj = UiC with C being 
the charge conjugation matrix. Since the Z-exchange amplitude vanishes at rest, only the latter 
amplitude contributes to the a coefficient of the resulting cross section (Jww^vei — Q>ww H~ ^ww^ e \- 



aww 



32tt 



1 



m 



w 



mi 



3/2 



1\ 



Vi x fml 



+ /i 2 -m 2 w ) 2 ' 



[mi 



(46) 



However, b-ww needs to be included since it turns out that in this case typically b ^> a and the 
correction to the cross section at freeze-out is of order 1. Both the Z-exchange and chargino- 
exchange amplitudes contribute to by/w- The resulting expression is very long, and we do not 
present it here. 

Finally, above the Z thresholds, we have to take into account the annihilation into ZZ pro- 
ceeding via i,u-channel exchanges of all the three neutralinos Xa {x = Xi) with the amplitude: 



2 3 
4cl V* 

a=l 



V2 rj ry 



m n 



(h -Pi) 2 -m 2 a 



7V«i -(1^2) 



e*(A;i)e*(A; 2 ), n a = V lx V 2a + V la V 2x . 

(47) 



The a coefficient of the resulting cross section azz can be given explicitly: 

-i 2 



azz 



9' 



647TC1 



m\ 



mi 



2 \ 3/2 



E 



2 

K a m x 



a=l X 



m 2 + m\ — m 2 z 



(4J 



However, bzz is again needed, which we do not present here. 

In our computation, we do not include the t-channel sfermion exchange contributions to the 
// final states or the s-channel Higgs-boson exchange contributions to the WW and ZZ final 
states, which are generically suppressed and depend on free parameters other than /i, M, A and 
tan/3. We also do not consider final states involving the Higgs boson(s), whose thresholds are 
higher than the ones discussed above. 
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tanp=1 .3 tanp=1 .4 




Figure 4: The lightest neutralino relic density Q x h 2 for four values of tan/3. The contours for 
Q x h 2 = 0.01,0.05 and 0.2 are denoted by the blue (thin solid) lines, while the blue (darkest) 
shading indicates regions with 0.09 ^ ^x^ 2 ~ 0.13, corresponding to the 95%CL region from 
WMAP [23]. Gray shading indicates regions where m x < mz/2. The contours for the mass of the 
LSP, m x = 80, 150 and 200 GeV, are also superimposed by dashed lines. 
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We now report the results of our numerical computations for the DM abundance. In Fig. 4 
we show plots for the DM relic abundance in the [x-M plane for 4 different values of tan (3 = t — 
1.3, 1.4, 1.5 and 1.6. In most of the parameter space for t pa (1 ~3), which we are interested in, 
the thermal relic abundance is much smaller than the observed DM abundance (< 10%). Around 
t = 1.4, there are regions where the abundance has the observed value. The plots show the 
situation for 4 representative values of t around 1.4. Going to smaller t is disfavored by the EWPT, 
while for larger t the DM region gets pushed towards larger values of /i, which are disfavored by 
naturalness. We conclude that the new possibility for DM discussed here is compatible with the 
EWPT and naturalness for tan/3 ^ 1.7. 

In the MSSM with decoupled gauginos, the lightest neutralino is a pseudo-Dirac Higgsino of 
mass \x. Such a Higgsino has a large coupling to the Z and, since naturalness places a bound 
on how large /i can be, the Higgsinos annihilate easily in the early universe contributing little to 
the DM. How does the LSP Higgsino of ASUSY avoid this? There are two possibilities. In the 
DM regions with low \M\, visible for the low values of tan/3 in Fig. 4, x is mainly singlet and 
this reduces the Z coupling. On the other hand, in the high \M\ DM regions, the Hi and H 2 
components of the Majorana x l ea d to a partial cancellation of the Z coupling, again allowing 
substantial DM. 

We finally consider the direct detection cross section of x- The x coupling to the lightest Higgs 
is given by 

xs 

~^XXh, 5 = y 3x (y lx cos(a-7r/4) + y 2x cos(a + 7r/4)), (49) 

while the coupling of h to the quarks is given by 

1 , / cos a _ sin a _ \ , . 

h ^—rm u q u q u ^m d q d q d . (50) 



s/2v V sm P COS P 

The nucleon cross section is then given by 

, . ... m 2 r ( \5 \ 2 ( cos a sinaA 2 2 
7r \m z h v J \ smp cosp/ 

where m r = m x m^ / (m x +mx) is the reduced mass, and X u and X d are certain linear combinations 
of nucleon matrix elements, which we conservatively take as X u ~ 0.14 and X d ~ 0.24 [24]. The 
cross section on a nucleus, normalized to a nucleon, is (for A = 2): 

/ \ aa 9 ^ ^ ,., / , r cosa ^ sin a A 2 / 300 GeVA 4 . , 

a h(xP ~^ XP) = 10 cm G, G~ 200 x§ 2 (X u ——-X d . 52 

\ sm jj cos [i J \ m h J 

Recall that for a DM mass in the 50 — 200 GeV range of interest to us, the limit on the spin- 
independent direct detection cross section is (2~3) x 10~ 43 cm 2 [25]. We have studied numerical 
values of the factor G in the region of tan/3 = 1.3-3 and found that it never exceeds ~ (20~30). 
Thus the cross section is not in contradiction with the existing limits, even without taking into 
account the fact that in most of the parameter space the LSP relic abundance is much below the 
observed DM density. We have further studied the factor G for the regions of parameter space with 
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tanp=1.3, m H +=400GeV tanp=1.3, m H ±=700 GeV 




Figure 5: Dark matter detection cross section a/, in units of 10~ 44 cm 2 (i.e. the factor G in the 
text) in the /i-M plane for tan/3 = 1.3, 1.5 and vnn± = 400, 700 GeV. The blue (darkest) regions 
indicate the regions for 0.09 ^ ^x^ 2 ^ 0.13 (see Fig. 4), while the gray shaded regions indicate 
the ones where m x < mz/2- 
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tan/3 = 1.3 - 1.6 where we get the observed relic density, i.e. the blue (darkest) regions in Fig. 4. 
We find that G ranges from about 1 to 5, as depicted in Fig. 5 for the case of = 400, 700 GeV 
for tan/5 = 1.3, 1.5. This puts the cross section (52) well within reach of experiments currently 
under way. 



5 Naturalness Constraints on Superparticle Masses 

As mentioned in the Introduction, making the lightest Higgs boson heavier helps in relaxing 
the naturalness constraints on the model, while allowing heavier superpartners. It is of interest 
therefore to establish on quantitative grounds the upper bounds that can be set on physical 
superparticle masses. This is of direct relevance also to understand the precise significance of the 
results obtained in the previous two sections. 

For this purpose, we consider an explicit simplified version of the ASUSY model, defined by 
the following superpotential and the soft supersymmetry-breaking Lagrangian 

W = XSHxHz + mH 1 H 2 + ^MS 2 , (53) 

Aoft = -m 2 s \S\ 2 - m iW 2 + (BmH 1 H 2 + h.c), (54) 

i=l,2 

and we will study especially the case in which the soft breaking mass mj is taken large, consistently 
with naturalness. (See Ref. [26] for an earlier analysis.) We do not expect the bounds obtained 
in the generic ASUSY model to be significantly different, or in any case more restrictive, than the 
ones valid in this specific case. 

The scalar potential takes the form of Eq. (4) with 

Hi = m 2 + /i 2 , h\ = m 2 + /i 2 , h\ = Bm — AMs, (55) 

V(S) = (M 2 + m 2 s )\S\ 2 =H 2 s\S\\ (56) 
and the chargino mass parameter defined in Eq. (3) is \x = m + As, where s is the background 
expectation value of S. In this model, the stability condition (6) implies also the conservation of 
CP at the vacuum. 

5.1 Minimizing the potential 

To understand the dependence of the minimum of the potential on various parameters, it is best 
to minimize first with respect to Vi and v 2 for fixed s. We then obtain 

V = His 2 - 1 (fas) - Hi(s)H2(s)) 2 . (57) 

This potential can in turn be minimized with respect to s, giving 

AV 

As = - 7 (/i 2 - H1H2) = - J - - —1 , (58) 



M + m(t + t~ 1 ] 
H 2 S + AV 



7= "r: v ;:: ; , m 
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Figure 6: The region with F > 1/2 in the \x-M plane, for t = 1.5 (between red solid lines) and for 
t = 3 (between blue dashed lines). 

where t = tan/3. Plugging this expression for s into Eq. (55) and then into (9), we can obtain 
the expression for v in terms of the original parameters of the model. Assuming that s is small 
relative to v, which is the case for sufficiently large fj,$, we find 

(4777, \ 
M + t + t _ 1 J 7, (60) 

which can be rewritten in the form 

AV ~ F' 1 [Bm (t + t' 1 ) - (m 2 + m 2 ) - (m 2 + m 2 )] , (61) 

(4m A 
M (62) 

This gives, together with t = fii/112, the explicit dependence of v on the original parameters of the 
model, which in turn allows us to determine the naturalness constraints. Note that the expression 
in the square bracket of Eq. (61) is the one which we would obtain in the pure 2HDM without 
the singlet field S. By taking m$ sufficiently large, consistently with naturalness (see below), and 
restricting the range of other parameters, we shall always require that the factor F be sufficiently 
close to unity so that we can neglect its presence. This is consistent with our treatment of the 
EWPT in Section 3, where the mixings between the singlet and doublet scalars were neglected. 
F > 1/2 is the numerical condition that we shall take to constrain the general space of the 
parameters. The regions F > 1/2 are depicted in Fig. 6 in the /i-M plane for t = 1.5 and 3. 
Here, the values of mn± and m| are taken to saturate the naturalness upper bounds discussed in 
the next subsection (see Eqs. (64, 66, 72)), although using smaller values of m H ± leads to only 
negligible variations of the regions. 
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5.2 Naturalness bounds 

We are interested in the naturalness constraints on /X3, which sets the bound on the Higgs sector 
through (11) and on the chargino mass through (20). We also want to know the bound on the 
stop masses and on the mass of the scalar S, which both affect, via one loop corrections, the 
masses m 2 and m 2 ., and therefore the masses fi 2 and /x 2 ,. As customary, we shall require that the 
logarithmic derivatives of v 2 with respect to fi 2 , fi 2 , mj~, and m| be less than A. 
One has, taking into account also the variation of tan j3, 

A V = 5/4 (t + r 1 ) - 5fil(2F! - F 2 ) - 5f4F 2 , (63) 

where 

(t - r 1 ) 2 / AV \ t 2 - 1 / \ 2 v 2 \ , . 

F 1 = l+ ± —L -z- + 1 , F 2 = 1 + — — -z- + 1 . (64) 



Here, we have used the approximation \i ~ m, /i 2 ~ Bm and F ~ 1, valid for small s. 

From the variation on /i 3 , upon use of Eq. (11), one obtains, taking A = 2 and normalizing to 
a fine-tuning A -1 = 20%, 

m A < 800 GeV(A/5) 1/2 . (65) 

This therefore sets, through (10, 20), the limits on the charged Higgs-boson and chargino masses 
(for A = 5): 4 

m H ± ^ 700 GeV, ji ^ cos (3 m H ±. (66) 

The stop masses, ttlq and m^, and the soft breaking mass for the scalar S, ms, affect mi and 
m 2 through the one loop RGEs 



dm 2 2 m 2 



+ (67) 



There is no particular problem in integrating these equations up to the messenger scale A. mess to 
obtain the standard stop contribution 

3 A 

6m 2 ~ \2^_2 , _2 M„ Yi mess 



»2 



l^ A 'K + m U ln TW' (69) 



with an RGE improvement needed if A mess gets far above 100 TeV or so. From the sensitivity of 
v 2 to /i 2 in Eq. (63), taking A mess = 100 TeV as a reference point, one obtains 



q £ 1.3 TeVsin/3(A/5) 1/2 /v/F 2 . (70) 



The contributions of m| to m 2 are, on the other hand, not equally well defined due to a rapid 
increase of A with energy. By integrating the RGEs up to A « 10 TeV, where perturbation theory 
is still valid, one obtains 

^(A^-^l + ^taA) „ -0.3ml, (71) 



4 The direct naturalness limit on fi from its contribution to ji\ and /i2 is less stringent. 
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where A represents the coupling A at A, and, consequently, 



m s < 1 TeV (A/5) 




(72) 



Finally, it is of interest to know the naturalness bound on the gluino mass rrig, which contributes 
to m\ via a two loop effect. Similarly to the bound on the stop masses, again for A mess = 100 TeV, 
one obtains 



Note that all these bounds are proportional to A, which is taken to be equal to 2. 

Figure 7 summarizes the knowledge of the spectrum for different values of tan f3 and for two 
reference values of m#± = 700,400 GeV. The former corresponds to the highest value compatible 
with the naturalness bound of Eq. (66), and both values are consistent with the constraint from 
b — > 57 without a destructive contribution from a stop-chargino loop [27]. The top and bottom 
panels show, respectively, the values of the Higgs boson masses and the upper bounds from 
naturalness on /i, m s , t R , and trig, allowed by a 20% tuning. The bound on the S'-fermion 
mass M is not explicitly given, since it depends, via F > 1/2, on many other parameters. 

We finally remark that the naturalness bounds derived from variations of tan j3 alone are less 
restrictive than the above bounds derived from the variations of v 2 . 

6 LHC Phenomenology 

How will LHC probe ASUSY — the theory defined by (3) and (4)? Clearly we must consider 
signals that result from the spectrum and mixings of the Higgs and Higgsinos, paying particular 
attention to the consequences of a large value of A, and to differences with the MSSM. A general 
and detailed treatment is far beyond the scope of this paper; we will be content with discussing a 
few illustrative examples, especially as they relate to a collider probe of DM. 

Throughout we assume that the S scalar is heavy and its mixing with the doublet scalars 
can be ignored. In this case, the parameters relevant for the doublet scalars are ii\ 2 3 and A, 
which are reduced to 3 by the minimization condition for the VEV v. Including both neutral and 
charged Higgsinos, in the limit that the gauginos are sufficiently heavy, there are an additional 2 
parameters fi and M, giving a total of 5 free parameters. It might appear that there are 2 extra 
parameters in this sector compared to the MSSM, A and M, but this is not the case — in the 
MSSM one needs an enhanced quartic coupling that brings in the top squark sector. Hence there 
is a single extra parameter, M, that describes the mass and mixings of an extra neutralino, so 
that the theory is highly constrained. Furthermore, the addition of the singlino allows for the 
possibility of Higgsino DM. 

The four doublet scalars themselves show important differences with the MSSM. First and 
foremost the lighter Higgs h is much heavier than in the MSSM. Second, the ordering of the 
spectrum is fixed: h, H , H with A heaviest. 5 In contrast, for the MSSM the ordering of H ± , 
A and H is not fixed, except A is lighter than H. Finally, tan j3 less than about 3 is excluded in 

5 Paramcter choices allow H ± to be lighter than h, but this must be consistent with the limit on m H ± from 
b — > S7. 




(73) 
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Figure 7: The spectrum of the Higgs bosons (top) and the naturalness bounds on parameters 
(bottom) as functions of tan/3, for two cases of m#± = 700 GeV (solid lines) and 400 GeV 
(dashed lines). In the top panel, four lines represent ttia, t^h, wih± and rrih from above, for each 
case of m H ± = 700,400 GeV. In the bottom panel, four lines represent m™ ax , m™ ax , rn™? and 
/i max from above, for each case of m H ± = 700, 400 GeV. 
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the MSSM, while it is strongly preferred in ASUSY, leading to important changes in the various 
decay rates for these scalars. The three neutral scalars h, H, A are each copiously produced at 
LHC from gluon-gluon fusion with cross sections in the range 0.1 - 10 pb, decreasing with mass. 
While these rates can be computed in terms of the Higgs masses and mixings, h, H and A may 
also be produced in cascade decays of top and bottom squarks (and from other squarks, if light 
enough to be produced) so that we will not consider the rates themselves for determining the 5 
free parameters, but will turn to the masses and widths of the three neutral scalars. 

With a branching ratio of about 10~ 3 in much of parameter space, h and H will be visible in 
"gold-plated" events where the scalar decays to ZZ, and each Z decays to 11, where I is an electron 
or muon. This will allow accurate measurements of and m#. Similarly, the decay A — > Zh will 
allow a measurement of tua from events of the form UUjj where each // reconstructs to a Z, as 
does the two jet system jj. 6 The measurements of rrih,H,A, together with the constraint from the 
electroweak VEV, will allow all 4 parameters of the Higgs sector to be determined. This includes 
both the quartic coupling A, so that one can evaluate the scale at which non-perturbative physics 
sets in, and tan/3 = Unlike the MSSM, the ratio of VEVs can be determined relatively 

easily. Notice that rrih,H,A, and hence the extraction of the 4 fundamental parameters, do not 
depend on any assumption about the superpartner spectrum. 

To extract /i and M, we specialize to the case that the gauginos are sufficiently heavier than 
the Higgsinos that we can ignore the Higgsino/gaugino mixing. In practice the accuracy with 
which the parameters can be extracted will depend on the size of the gaugino components of 
the lighter neutralinos and charginos. Over much of the parameter space, the decay of H, A 
and h to Higgsino pairs is kinematically allowed. Since it is induced at tree-level by the large 
coupling A, it will be competitive with the WW mode, generically having a significant effect 
on the total widths Th,a,h- Note that the WW, ZZ and tt modes depend on the Higgs mixing 
parameters and on tan/3 (needed for example to compute the top Yukawa coupling). However, 
they can be accurately computed, since these parameters will have been determined by the mass 
measurements. The widths to Higgsinos depend on both the Higgs mixings, already determined 
by the Higgs masses, and by the Higgsino mixings and masses, which are determined by \x and M. 
Can the total widths Tn,A,h be measured sufficiently accurately, by reconstructing the invariant 
mass from gold-plated events, to extract fi and Ml For A — > Zh there will be insufficient rate to 
rely on the 61 events, and the 41 jj events are unlikely to be reconstructed with sufficient accuracy. 
Hence we consider the possibility of measuring \x and M via the contributions of the Higgsino pair 
mode to IV and I\, reconstructing these total widths using the 41 final states. 

The partial widths to Higgsinos for H and h are shown as contours in Fig. 8 for A = 2, 
tan (3 — 1.3, 1.5 and m#± = 400, 700 GeV in the /i-M plane. For each of the four panels we give the 
dominant non-Higgsino widths for comparison: T (h) = T(h — > WW + ZZ) and T (H) = T(H — > 
WW + ZZ + tt + hh). Consider for example the panel with tan/3 = 1.3 and m#± = 400 GeV, 
where To(h) ~ 8 GeV and To(H) « 6 GeV. From inspection of the contours, we find that the 

6 As m,H increases and approaches its naturalness limit of 700 GeV the decoupling regime is reached, so that 
the decays H — > ZZ and A — > Zh become suppressed relative to the tt mode. Since the production cross section is 
also reduced, mass measurements from the gold-plated decays may be difficult in this region. However, this region 
is not the most natural expectation. Note that while the couplings of either h or H to vector boson pairs vanish 
in ASUSY as tan/3 — > f, this suppression becomes rapidly ineffective as tan/3 starts deviating from f. 
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Figure 8: Contours for the partial widths H — > Higgsinos (dashed red curves) and h — > Higgsinos 
(solid black curves), labeled in GeV, in the fi-M plane for A = 2. The four panels are for 
tan/3 = 1.3,1.5 and = 400,700 GeV. For each of the four panels we give the dominant 

non-Higgsino widths, To, for comparison. The gray shaded regions have m x < m^/2, while the 
dark (blue) shaded region corresponds to 0.09 ^ Q x h 2 ^ 0.13. 
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Higgsino widths for both h and H are in the range of (10 ~ 100)% of the corresponding non- 
Higgsino widths, allowing a determination of /i and M over much of the parameter space. The 
accuracy of the determination will depend on the location in the [i-M plane. For much of the DM 
region, shown shaded in Fig. 8, the percentage changes in the widths induced by the Higgsino 
mode is roughly in the region of 30%. Increasing m H ± to 700 GeV leads to similar percentage 
changes in the widths, so again measurements of fi and M are possible. The number of gold-plated 
events from H decay decreases as both rrtu and tan j3 are increased, making the measurement more 
difficult as the regions disfavored by naturalness and EWPT are approached. It will, of course, 
be highly significant if measurements at the LHC indicate values of \x and M in the DM region, 
determining whether the DM Higgsino is mainly singlino or has comparable doublet and singlet 
components, and allowing predictions for the direct detection rate and the LSP mass. As tan (3 is 
increased from 1.3 to 1.5 the signals in Th,u are comparable, but the dominantly singlino region is 
essentially absent, so it will be an easier task to verify consistency with the remaining DM region. 

In Fig. 9 we show the contours for T(H, h — > Higgsinos) for larger values of tan/3 that do not 
allow thermal relic Higgsinos to be DM. Nevertheless, the change in the total widths induced by 
the Higgsino contributions is again a significant percentage, so that [i and M can be measured. 

Determining parameters from studying the decay modes of Higgsinos may also be possible, 
but it looks to be very difficult. Consider pair production of top squarks, followed by decays to 
t + Higgsino. The heavier neutral Higgsinos cascade to the LSP and a Higgs boson, so that the 
decay chain is t — > t\2,z followed by x 2 ,3 — ► XiH or \2,z XiA leading to signals tillll with the 
four leptons reconstructing to or rrif L . The ratio of these two types of events is another probe 
of the parameter space, since it is independent of the mass and production rate of the top squarks. 
However, this could be contaminated by events with bottom squarks decaying to top quarks and 
the charged Higgsino, which itself cascades to H or h. 

7 Overview and Conclusions 

One year from now the Large Hadron Collider will start directly exploring for the first time 
the energy range well above the Fermi scale. Many speculations of the last three decades on 
Electro Weak Symmetry Breaking (EWSB) will be rendered irrelevant, while some, perhaps, may 
emerge as part of physical reality. Therefore, the time is right to briefly reconsider the status of 
these speculations, including the one presented in this work. 

There are two main sets of considerations that have played a dominant role to orient ideas 
and specific work on the subject. One rotates around the naturalness problem of the Fermi scale. 
Why is the Higgs boson light relative to any mass scale appearing in the UV completion of the 
SM, whatever it may be? More quantitatively, as illustrated by all known examples that attack 
the problem, why is the Higgs mass not at least as large as the SM contribution with a cutoff 
scale where new physics sets in, Anp? The second set of considerations relates to the significance 
of the EWPT to the EWSB problem. If new physics occurs at Anp, would it not be manifest, 
even indirectly, in EWPT? 

In principle these two kinds of considerations do not conflict with each other: they might have 
actually merged into a coherent picture for the physics that underlies EWSB. So far, it is fair to 
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Figure 9: As for Fig. 8, except for tan/3 = 2,2.5. There are no regions giving sufficient thermal 
relic Higgsino DM. Note that a larger region of M is shown than in Fig. 8, but smaller regions in 
fi, corresponding to the more stringent naturalness constraints at larger tan (3. 
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say that this has not happened. On the contrary, the impressive success of the SM in providing 
an accurate perturbative description of the EWPT casts doubts on the physical relevance of many 
ideas put forward to attack the naturalness problem of the Fermi scale. This issue, the "little 
hierarchy problem," is a quantitative one [28]: for A^p as low as implied in various attempts 
to attack the naturalness problem, why is the new physics not manifest in the EWPT? At first 
sight the supersymmetric attempts seem exempt from this difficulty, as they cope easily with the 
EWPT. In their simplest versions, however, they also have a problem accommodating the direct 
bound on the Higgs mass. To satisfy the bound apparently requires either a lack of naturalness 
or a complication of the theory. 

It would neither be possible nor appropriate to recall here all the attempts to get around this 
difficulty. Nevertheless one aspect is worth emphasizing, as trivial as it may be. While ignoring the 
naturalness problem is not in contradiction with experiment, it has however a striking practical 
consequence, other than being theoretically unsatisfactory: it weakens the case for finding new 
physics at the LHC. Needless to say, what will be found at the LHC does not depend on our choice 
of the physically relevant problems. However, insisting on the naturalness problem, as we do in 
this paper, has at least the advantage that the proposals to solve it will be scrutinized when the 
LHC will be successfully turned on. 

With this motivation in mind, the framework analyzed in this paper originates from the obser- 
vation that the little hierarchy problem might have (part of) its source in a misinterpretation of 
the EWPT, especially related to the inferred strong bound on the Higgs mass. Maybe the results 
of the EWPT, although consistent with the SM, do hide after all some new physics. This might 
be especially the case if the corresponding new theory still allows a successful perturbative de- 
scription of the EWPT in a large portion of its parameter space. At first sight, this does not seem 
to be of relevance to supersymmetry, since, as observed above, standard supersymmetry has no 
special problem with the EWPT. However supersymmetry in its minimal version has a theoretical 
bound on the Higgs mass which is even stronger than the indirect experimental bound from the 
EWPT, causing its own naturalness problems. Hence the proposal that we make here of relaxing 
as much as possible the upper bound on the Higgs mass by going to ASUSY, while still retaining 
the successful perturbative description of the EWPT. In our view, this "bottom-up" approach to 
the little hierarchy problem deserves attention. 

We can summarize this bottom-up approach to supersymmetry as follows: the Higgs sector is 
pushed up to the (200 - 700) GeV range, the superpartners may then be heavier than 700 GeV, 
while most of the parameter space with successful EWSB naturally yields the observed values for 
the S and T parameters. 

Needless to say ASUSY has two related prices to pay, which have certainly restrained most 
people from considering it seriously up to now. The first is that it is not a UV-complete theory in 
the same sense that the MSSM or the standard NMSSM are: at an energy scale above 10 TeV or so, 
some change of regime of the theory should intervene that prevents the coupling A from exploding. 
This is not uncommon in most attempts to address the little hierarchy problem. The other price is 
that ASUSY is not manifestly consistent with gauge coupling unification. To determine this would 
require specifying the form of the necessary change of regime at high energy; it is nevertheless 
remarkable that a successful solution already exists in the literature [10, 11, 12]. 
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A One-loop Contributions of New Particles to S and T 

Recall that the T and S parameters are given by 

A 33 (0) - A ww {0) 
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where Am and Fij are extracted from the gauge-boson pair vacuum polarization amplitudes: 



nJU?) = -ig^iO) + q 2 F l >(q 2 )} + g^terms, 



(76) 
sin 9 W 



with i,j = W, 7, Z, or i,j = 0, 3 for and W^^, a cm is the structure constant, and s w 
and c w = cos9w with 0w the Weinberg angle. 

We record here one- loop expressions for A(0) [29] and F(0) produced by bosons and fermions 
coupled to a generic gauge boson with unit strength. For a boson loop with internal masses 

mi and m 2 and coupling iW^l d M 02 we have 
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There is an additional diagram contributing to S for a Higgs boson, in which the gauge boson 
and the Higgs boson <p propagate in the loop, giving SS = m^Gim^^mw) where 
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For a fermion loop with internal masses mi and m 2 and a vector coupling Wu&i'Y^fai we have 
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For an axial coupling, the result can be obtained by letting m\ — > —mi. These expressions are 
valid for both Dirac and Majorana fermions, with an extra factor of 2 in the case of identical 
Majorana fermions. 
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